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Abstract. We present strong approximations with rate of convergence for 
the solution of a stochastic differential equation of the form 

dXt = b(Xt)dt + (T(Xt)dBf , 

where b G C^, o G C^, is fractional Brownian motion with Hurst in- 
dex i?, and we assume existence of a unique solution with Doss-Sussmann 
representation. The results are based on a strong approximation of by 
means of transport processes of Garzon et al (2009). If a is bounded away 
from 0, an approximation is obtained by a general Lipschitz dependence re- 
sult of Romisch and Wakolbinger (1985). Without that assumption on ct, 
that method does not work, and we proceed by means of Euler schemes on 
the Doss-Sussmann representation to obtain another approximation, whose 
proof is the bulk of the paper. 



1. Introduction 

We consider one dimensional fractional stochastic differential equations of the 
form 

dXt = b{Xt)dt + a{Xt)dB^ , iG(0,T], (1.1) 
Xo = Xo, 

where is fractional Brownian motion with Hurst index H, and b and a are 
continuous functions. Equations of this type appear in several areas of application 
due to the properties of (see e.g. [18] and references therein). B^ = {B^)t>o 
is defined for any H G (0, 1) as a centered Gaussian process with covariance 

Co.{B^,B^)=^-is'"+t'"-\s-tn. 

We exclude the case H — 1/2, which corresponds to Brownian motion. The main 
properties of B^ (for H ^ 1/2) are self-similarity, stationarity of increments, 
long range dependence, fc-H61der continuity of trajectories for k < H, and it is 
neither a Markov process nor a semimartingale, hence the classical Ito calculus 
cannot be used for this process. See [23, 27] for background. There is now an 
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extensive literature on fractional stochastic differential equations and applications. 
We restrict the references to a minimum for reasons of space. 

In order to give a precise meaning to the stochastic differential equation (1.1) it 
is necessary to define stochastic integral with respect to and to define solution, 
and then existence and uniqueness of solution can be proved under suitable condi- 
tions on H, b and a, in particular for the solution to have a Doss-Sussmann repre- 
sentation (given by equation (2.9)). For example, this can be done for H > 1/3, b 
Lipschitz, and a G C^. These questions have been treated in [21] (see also [4] and 
references therein). The paths of are increasingly irregular as H decreases, 
and H = 1/4 is a critical value for some problems (see e.g. [5, 10, 14, 15, 22]). 
Hence the assumption iJ > 1/4 is often made. In this case, under some conditions 
on b and cr, a unique solution with Doss-Sussmann representation has been estab- 
lished in [1] for 1/4 < < 1/2 with Stratonovich integral, and in [20] for H > 1/2 
with forward integral. In [2] it is remarked that under some assumptions on the 
integrand, the Stratonovich integral coincides with the forward integral (see [26] 
for those integrals). In this paper we take b G and a G C"^ {CI is the space 
of bounded functions with continuous bounded derivatives of order < i), and we 
assume existence of a unique solution of (1.1) with Doss-Sussmann representation; 
this holds at least for H > 1/4. The restriction i? > 1/4 is not used in our proofs, 
and it seems to be unnecessary if (1.1) is treated with the theory of rough paths; 
moreover, in that way it may also be possible to study strong transport approx- 
imations in the multidimensional case (see [3, 17, 24, 28] and references therein 
concerning the theory of rough paths and related fractional stochastic diferential 
equations) . 

In [11] we obtained a strong approximation of B^ with a rate of convergence by 
means of the Mandelbrot-van Ness representation of B^ as a stochastic integral 
with respect to Brownian motion [16], and a strong approximation of Brownian 
motion by transport processes [13]. The result is that for each q > and each (3 
such that \H — 1/2| < P < 1/2, if is the n-th approximation of B^ (defined in 
(2.6)-(2.7)), there is a positive constant C such that 



Thus, the approximation becomes better when H approaches 1/2. 

The aim of this paper is to show that under the above assumptions on b 
and CT, approximate solutions X" of equation (1.1) can be constructed by means 
of S", which strongly approximate the solution X similarly as (1.2), with rate 
n~-'^/^"'"''"'"''(logn)^/^, where (5 > is arbitrarily small, i.e., almost the same rate 
as (1.2). If <j is bounded away from 0, this can be done directly from (1.2) us- 
ing a general Lipschitz dependence result obtained in [25] (see also [8]), without 
recourse to specifics of B". This is Theorem 3.1. Without that assumption on <j 
the result of [25] cannot be used (see Remark 4.1), and we proceed by applying 
Euler schemes to the Doss-Sussmann representation of the solution. In this case, 
properties of _B" arc involved, and the approximation is somewhat different in the 
probability sense. This is Theorem 3.3, whose proof is the bulk of the paper. 




as n — > oo. 



(1.2) 
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We stress that the approximations of the solution of (1.1) by means of trans- 
port processes are of theoretical interest, but may not be useful for generating 
approximate solutions with computational efficiency. To that end there are ad hoc 
methods (e.g. [19] and references therein). 

Section 2 contains background, Section 3 results, and Section 4 proofs. 



2. Background 

2.1. Approximation of fractional Brownian motion by transport 
processes. 

For each n = 1, 2, . . ., a (uniform) transport process = {Z'^{t))i>o represents 
the position on the real line at each time t of a particle that starts from with 
velocity +n or — n, with probability 1 /2 each, continues with that velocity during 
an exponentially distributed time with parameter n^, at the end of which it changes 
velocity from ±n to =Fn, and so on the same way, changing sign at consecutive 
independent exponentially holding times. Such a process can be constructed from 
a given Brownian motion i? on a probability space using the Skorohod embedding, 
and it was shown in [13] that converges strongly to B uniformly on a given 
bounded time interval with rate Cn~^/^ (log n)^/^, as n — > oo, where C is a positive 
constant. 

The Mandelbrot-van Ness representation of is given by 

B^ = Ch(^J° ft{s)dB{s) + j\t{s)dB{s)y 0<t<T, (2.1) 

where 

f,{s) = (t - s)^-i/2 _ (_s)//-i/2 for s<0<t, (2.2) 

gt{s) = (t-s)"-'/^ for s<t, (2.3) 

B is a Brownian motion on the whole real line, and Ch is a positive constant [16]. 
Fix a < 0. After an integration by parts and a change of variable, B^ can be 
written as 

= Ch ( f 9t{s)dB,{s) + r Ms)dB2{s) + Ma)B2{a) 

\J0 Ja 

djt (-) ^B3{v)dv] , (2.4) 



V J 



n/a 

where Bi , B2 and B^ are Brownian motions given, respectively, by the restriction 
of B to [0, T], the restriction of B to [a, 0], and 



sB(i) if sG[l/a,0), 
if s = 0. 



To define the approximation S" = (B")t>o that appears in (1.2), the idea is to 

approximate Bi,B2 and B^ by corresponding transport processes Z", Z2 and Z^. 
Forn = 1, 2, . . . and < ^ < 1/2, let 

£„ = -n-'5/l«-i/2|. (2.5) 
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For if > 1/2, define 

= Ch( f 9tis)dZ^{s) + f ft{s)dZ^{s) + ft{a)Z^{a) 



'1/ 

and for < 1/2, define 



/ gt{s)dZns)+ / gt{en + s)dZ^{s) 

Jo J(t+£„)VO 

+ Ms)dZ^{s) + Ma)Z^{a) + £ djt ^d^;) dZ^{s)j . 

(2.7) 

Note that depends on through (2.5). 

The following consequence of (1.2) is obvious since q is arbitrary, hence large 
enough. 

P (limsup I sup |Bf -B^\> Cn-i/2+''(logn)5/4 ) = 0' (2-8) 

V n^oo {o<t<T J / 

where lim sup is understood in the sense of sets. This will be used for the proof of 
Theorem 3.3. 

2.2. Doss-Sussmann representation. 

Under suitable assumptions on b, a and H, the Doss-Sussmann representation 
of (1.1) is given by (see [1, 7, 20, 21]) 

X, = h{Yt,B^), (2.9) 

where the function h and the process Y are the solutions of equations 
dh 

{X1,X2) = (T{h{xi,X2)), h{xi,0) = Xi, Xi,X2€M., (2.10) 

OX2 

and 

= exp(^- J^' a'{h{Yt,s))ds^b{h{Yt,Bi')), Yo = xo, (2.11) 
respectively. The function h has the property (see [7], Lemma 2) 



dh 

which implies that 

dxi 



^^^{xi,X2)=exp( I a'{h{xi,s))ds] , (2.12) 



y/= (|^(r„i?f)) b{h{YuB^)), Yo = xo. (2.13) 



APPROXIMATIONS OF FRACTIONAL STOCHASTIC DIFFERENTIAL EQUATIONS 5 



3. Results 

Recall that b e C^, a e C'^, (2.9) holds, and is defined by (2.6)-(2.7). 

The first theorem is a special case with the assumption that the function a is 
bounded away from 0. The interest of this result is that a direct proof can be 
given using Fernique's theorem [9], and a general Lipschitz dependence result of 
[25], without involving anything special about i?". 

Let X"- be the solution of (1.1) with replaced by S", and the integral is 
defined pathwise. 

Theorem 3.1. Assume a is bounded away from 0. Let — 1/2| < (3 < 1/2. Then 
for each 6 > such that /3 + 6 < 1/2, there exist q > and a positive constant C 
such that 

p( sup \Xt-X^\>Cn-^/'^+'^+^{\ognf/A =o{n-'') as n ^ oo. (3.1) 

Remark 3.2. (1) A transport approximation for equation (1.1) with Brownian 
motion {H = 1/2) and a bounded away from was studied in [12, 25]. 
The result is like (3.1) with /3 = 5 = 0. There was an error in the proof in 
[12], which was remedied with the method of [25]. The paper [6] gives a 
related result with a diflFerent type of formulation. 
(2) The approximation (1.2) holds for arbitrary q, but in the approximation 
(3.1) 5 is arbitrary and q is chosen appropriately small (see the proof). 
Hence a result like (3.2) below may not hold under the assumption of 
Theorem 3.1. 

For the general result, let X" — /i"(y"." ^ B"), where (/i")n is an Euler scheme 
approximation of equation (2.10), and (F" ™)™ is an Euler scheme approximation 
of y", which is the (Doss-Sussmann) solution of equation (2.11) with replaced 
by B^ (the precise definitions of ft," and y"-™ are given in (4.3) and (4.9)). 

Theorem 3.3. Let \H - 1/2| < /3 < 1/2. Then for each 6 > such that 5 < fi 
and /3 + (5 < 1/2, there exists a positive constant C such that 

P ( limsupj sup \X, - Xl'\ > Cn-^/'^+l^+\\ognY'/H) = 0. (3.2) 

V n^oo [o<t<T J / 

Note that without the assumption on a in Theorem 3.1 the same rate of con- 
vergence holds, but the probabilistic part of the result is slightly different. 

4. Proofs 

Wc write II • I loo for the sup norm on [0,T]. 

4.1. Proof of Theorem 3.1. 

We give an idea of the proof. By Fernique's theorem [9] there exists a > such 
that 

P{\\B"\\^>s)<Ke-'''\ s>Q, 
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where K = Eex.'p{a\\B^ W^) < oo. Then, applying Theorem 3 A) and Theorem 4 
of [25], from (1.2) we obtain 

P{\\X - > n'^ Cn-^/2+^ (log n)5/2) = o(n-«) as n oo, 

with S = qj/a^^^, whore 7 is the positive constant in Theorem 4 of [25]. Since q 
in (1.2) is arbitrary, it can be chosen so that 5 is as small as desired, and (3.1) is 
obtained. □ 

Remark 4.1. It can be shown that for (t(.t) = arctan x, formula (40) of [25] cannot 
be proved because 

— (x,2/) >exp' 



dx ^ J \ 1 + x' 

for a; < 0, y < 0, so (§|(a:, y))~^ is unbounded, and therefore the argument for the 
proof of Theorem 4 in [25] is not valid. 

4.2. Proof of Theorem 3.3. 

Since the constant C in (2.8) does not play an essential role, for simplicity we 
put C = 1, and we will also prove the theorem with C = 1 in (3.2). 

Wc start by describing the parts of the approximation. We first approximate 
the function h and the process Y given by (2.10) and (2.11) respectively, and then, 
based on those approximations, we formulate the approximation for the solution 
of (1.1). 

The function /i : _^ ^ satisfies 

i-y 

h{x,y)=x+ / a{h{x,s))ds. (4.1) 
Jo 

For each n — 1, 2, • • • , wc take the partition {j/"} of the interval [— n, n] given by 
-'^ = V-n^ < • • • < < ?/o = < 2/1" < • ■ ■ < 2/^2 = n, where for r„ = 1/n, and 
i = 1, • • • , — 1, 

yr+i=i/r + r„ = ^, yVi)=y-i-^" = -^- (4-2) 

We define the functions ft." : _^ M by 

h"'{x,y) = if (.T,y) ^ [-n,n] X [-n, n], 

for {x, y) G [— n, n] x [— n, n] and A: = 0, 1, • • • , — 1, 

h-{x,y^) = X, 
K^{x,yl+^) = h^{x,y^) + r^aih"ix,y^)), 

h^{x,yl^k+i)) = h^{x,ylk)-rn(T{h^{x,ylk)l 
and by linear interpolation, 
h^{x,y) = h^{x,yl) + {y-y^)a{h^{x,yl)) if yl<y<yl+^, 
h^{x,y) = h^{x,yl^)-{ylk-y)cj{h^{xM) if yl^^+i) < V < V-k- (4-3) 
From (2.13) we have 

Yt = xo + £ (^^(n, i3f )) 6(/i(n, B^))ds. (4.4) 



APPROXIMATIONS OF FRACTIONAL STOCHASTIC DIFFERENTIAL EQUATIONS 7 



For each n = 1, 2, • • • , we define the process Y" as the solution of 

Yr = xo+J^(^§^^{Y:,B:)^ \{h{Y:,B^))ds, te^. (4.5) 

These processes exist and are unique because for each n= 1, 2, • • • , the function 
{x,t)^(^^^{x,B^)^ b{h{x,Bt)), t€[0,T], 

satisfies Lipschitz and Unear growth conditions in x (sec (2.12), (4.7) and Corollary 
4.5 below). Moreover, they are defined on the same probability space as the 
Brownian motion B in (2.1), since the solution of (4.5) is given for each sample 
point, and B" is defined on the same probability space as B. 

Wc give next an Euler scheme for approximating the processes F" for each 
n=l,2,---. 

By (2.12), equation (4.5) can be written as 

{Yn' = f{Y-,B-), Y- = xo, (4.6) 
where f{x,y) is defined by 

f{x, y) = exp a'{h{x, u))du^ b{h{x, y)). (4.7) 

For each n = 1, 2, • • • , we define /"(x, y) by 

r{x, y) = exp cj\hr{x, u))di^ h{hr{x, y)), (4.8) 

with /i" given by (4.3). 

The Euler scheme (F"'™)™ for equation (4.6) is defined as follows for each 
m = 1, 2, • • • , the partition = to < ■ ■ ■ < tm = T oi [0,T] with ti+i = U + rm, 
and rm — T/m: 



y'Ti.m 

=2^0, 

YtlZ = Y^r+rmFiYtT^K)' fc = 0,...,(m-l), 

yn,m ^ Y,-^^ + {t-tu)r{Yrr,Bl) 



(4.9) 



Equation (1.1) has a unique solution X with representation (2.9), and wc define 
similarly the approximation X" by means of the Euler schemes for h and as 

x;^ = /i"(y,"'"',Bn> (4-10) 

where /i" and F"'"^ are given by (4.3) and (4.9), respectively. 
We then have 

\Xt - I < H^{t) + H2{t) + H^{t), (4.11) 
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where 



\hiYuB^)-hiYr,B-)l 
\h{Yr,B-)-hiYr\B-)\, 

\h{Yr\B^)-h-{Yr\B2)\. 



(4.12) 
(4.13) 
(4.14) 



The proof consists in obtaining estimates involving Hi , H2 and H3 . 

We will need the following preliminary results. First a Lipschitz property of 

Lemma 4.2. Let B"- be defined by (2.6)-(2.7). Then for all n and for ti,t2 G 
[0,T], 

\Bl-Bl\<Kn^+^ \t2-ti\, 
where K is a positive constant. 

Proof. First we take H > 1/2. It suffices to check the property for each one of the 
four functions on the r.h.s. of (2.6). using (2.2), (2.3) and the transport processes. 
We will omit some of the calciilations. 

Note that dgft > by a straighforward calculation (see [11], Lemma 3.1). 

For ti <t2, by the mean value theorem, 

ft2 



Jo Jo 

(t2-sf-'/'-{ti-sf-'/' 



ds 


+ 


J 






Ju 



\t2-sf-'/'ds 



H+1/2 

Again by the mean value theorem, 
,.0 rO 



[ ft,{s)dZ^is) - [ ftMdZq{s) 
<n 1° 1(^2 - s)"-^/^ - {ti - s)"-y^\ds 



<- 



■H+1/2 



{{t2 - af+'/^ - (ti - a)^+i/2) < n{T - af-''^ \t2 - ti| . 



Next, 



\ftM)Z^{o.) - /ti(a)^2 («)l < n{-a)\{t2 - a)"-''^ - (ii - a)^-!/^! 
<n{-a)"-^'^{H -1/2) \t2-ti\. 
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Finally, 



pO / psAe 
Jl/a\ Jl/a 



i\ 1 



dv]dZ!i(s) 



'l/o\ Jl/a J Jl/a\ Jl/a 

< n{H - 1/2) / / I - (^2 - + (ti - l/v)"-^/^\{-l/vfdvds 

Jl/a Jl/a 

< n{H -l/2){3/2- H)\t2-ti\ / {-v)-"-'^/'^dvds, 

Jl/a Jl/a 



and, by (2.5), 



Jl/a Jl/a n - l/Z Jl/a 



1 



if - 1/2 



hence the result for H > 1/2 follows. 
We proceed similarly for H < 1/2. Let 



L 



(t2+£n)V0 



9t2 is)dZ- 



lis)- I 
Jo 



(tl+£n)VO 



If + £„ < and t2 + Sn < 0, then Ai =0. If + £„ < and ^2 + > 0, then 
by the mean value theorem. 



Ai 



{t2 - s)"-'/^dZ^{s) 



<n {t2- s)"-^'^ds < n{-enf-'/^{t2 + Sn) < n'^^ \t2 - h\ 

Jo 

and if ti + £„ > and ^2 + £n > 0, again by the mean value theorem. 



Ai < n 



tl+En /■*2+eTi 

Jti+e„ 



< 2n(-£„)^-i/2 \t2-h\= 2n'+^ \t2-h\. 



Let 



/ 5*2 (£n + s)dZ'l{s) - / gt^ (£„ + 

•'(t2+e„)V0 i(ti+£„)VO 



■s)(iZi"(s) 
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We show only one case. The others are similar. If ti + Sn < 0, t2 + Sn > Q and 

tl >t2 + £n, 



tl 



't2+S 



t2+e„ 



< 



2n 



(tl-£„-s)«-l/2dZi"(s) 
r [{t2 -Sn- s)"-y^ - (ii - e„ - s)"-y^]dZns) 



+1/2' 

< 2n(-£„)^- -ti\ = 2n'+f^ \t2 -h\. 



Next, 



<- 



J a 



f 

J a 



{t2-en)"+'/'-{h-er.)"+'/' 



■H + 1/2 

<n{h - £n)^-'/' \t2 - til < n(-e„)«-i/2 \t2 - h\ = n'+^ \t2 - h\ , 

IftMZ^a) - ftM)ZnO')\ < ni-a){l/2-H){h-a)"-^/'\t2-h\ 

< n'+^i-a)"-'/\l/2-H)\t2-h\, 

I-(/>'^''(^)^''")^^?'''-/l(-/>-^*-(^)>')^?<'' 



<n 



l/a Jl/a 



{1/2 -H) ih-l/v) 



l/a \ Jl/a 



<n / (1/2 - H){3/2 - H){ti - l/vf'^/^ \t2 - h\ {-1/vfdvds 

Jl/a Jl/a 

<(3/2 - H){-l/af/^-"n^+^ \t2 -ti\. 



□ 



Corollary 4.3. For each n = 1, 2, • ■ • , < oo a.s. 

Proof, immediate from Lemma 4.2. □ 
By the assumptions on a and b, we have the next bounds: 

\b{x)\ < Ml, \<j'{x)\ < M2, \a"{x)\ < M3, 

- b{y)\ < M4,\x - y\, \b'{x)\ < and |cr(a;)| < M5, (4.15) 

where Mi, • • • , M5 are constants. We put M = max{M2, M5}. 
Next we present some properties of the function h. 

Lemma 4.4. Let h be defined by (2.10). Then 
(1) 



dh 
dxi 



{x,y) 



< exp(M2 \y\). 



(4.16) 
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(2) 



(3) 



(4) 



exp 



fV 

Jo 



X, s))ds 



< exp(M2 \y\). (4.17) 



d ( dh 



<M3|y|exp(2M2|y|). 



(5) 



dxi 



dxi \dxi 

\h{xi,y) - h{x2,y)\ < exp(M2 \y\) \xi - X2 
dh 



5r-(a^i'y)) -(£^(^2,2/) 



9a; 1 



(4.18) 
(4.19) 

<M3|y|exp(2M2|y|)|a;i -a;2|. (4.20) 



(6) 
(7) 
(8) 



\bihixi,y)) -bih{x2,y))\ < M4exp(M2 |y|) -X2I 
\b {h{x, yi))-b {h{x, 2/2)) | < M4M5 |yi - ^2! • 



(4.21) 
(4.22) 



dh 
dxi 



{xi,y)) b{h{xi,y)) 



dh 

dxi 



ix2,y)] b{h{x2,y)) 



<exp(2M2 |t/|)[MiM3 \y\ + M4] [xi - X2I . 

Proof. (1) and (2) arc immediate from (2.12) and (4.15). 
(3) As a" is bounded by Afg, and by (4.16), (4.17), then 



(4.23) 



d f dh 



dxi \dxi 



■ix,y) 



exp ^— y a'{h{x,u))du 

< M3|y|cxp(2M2|y|). 



y gfi 
a"(h(x, u))— — (x, u]du 
0x1 



(4) Follows by the mean value theorem and (4.16). 

(5) Follows by the mean value theorem and (4.18). 

(6) As b is Lipschitz, the result follows by the mean value theorem and (4.19). 

(7) Using (2.10), b being Lipschitz and a bounded, for some ^ between yi and 
2/2 we have 



dh 

dX2 



\b{h{x,yi))-b{h{x,y2))\ <MA\h{x,yi)-h{x,y2)\=Mi 

= M4a{h{x,^))\\yi-y2\ < M^M^lyi - y2\. 

(8) 



\yi - 2/2 1 



dh 



dxi 



— {xi,y)] b{h{xi,y))- [—{x2,y)] b{h{x2,y)) 



dh 



dxi 



<Ai+A2, (4.24) 
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where 



dh 
dxi 



dh 
dxi 



\h{h{x^,y))\. 



A2 = \b{h{xuy))-b{h{x2,y))\ 



dh 
dxi 



{x2,y) 



Using h bounded by Mi, (4.21), (4.17), (4.20) and (4.24) we obtain the 
result. 

□ 

Corollary 4.5. Let f be defined by (4-7). Then for each n = 1, 2, • • • the function 

{x,t)^f{x,Bl^) 
satisfies Lipschitz and linear growth conditions in x. 

Proof. By (2.12), Lemma 4.2 and parts 1 and 8 of Lemma 4.4 , for each xi,X2 € M 
and t e [0,r], 

\f{xuB^) - f{x2, B^)\ < exp{2M2Kn^+^T)[MiMsKn^+'^T + M4] \xi -X2\, 
and as b is bounded by Mi, 

\f{xuB^)\ < MieMM2Kn^+^T). 



□ 



Lemma 4.6. Let f be defined by (4-7). Then for each n = 1, 2, • 

(1) |/(a;,i?r)|<MiCxp(M2||i?"||^), 

(2) \fixi,B^)-fix2,B^)\<Zi\xi-X2\, 

(3) \fix, Bl) - fix, BZ)\ < Z2 \B^, - Bl\, 
with the random variables 

Z, = exp(2M2 ||i?"|L)[MiM3 + M4] 

Z2 = (M1M2 + h'hMi) cxp(M2 ||B"IL)- 

Proof. (1) Follows from (2.12), (4.17) and boundedness of b. 

(2) Follows from (4.23). 

(3) By (2.10), 

df{x,B^) 



(4.25) 
(4.26) 



dX2 



f{x,B^)a'{h{x,B^)) 



+ exp 



a'{h{x,u))du b'{h{x,B^))(7{h{x,B^)). 



Since a, a' and b' are bounded, and by part 1 and (4.17), then 
df{x,B:^) 



8X2 



< (M1M2 + M5M4) exp(M2 



and we have the result by the mean value theorem. 



□ 
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Now we do the approximation of h. For fixed n, we work in the square [— n, n] x 
[— n, n]. Let I = n + m for some m > 0, and consider the finer partition of [— n, n] 
given by — n = j/L„; < ■ • • < 2/o = < • • • < j/^; = n, as in (4.2) with r; = l/l. 

Lemma 4.7. Let h and be given by (4-1) o-nd (4-3), respectively. Then for 

{x, y) € [— n, n] x [— n, n] and I > n, 

\h{x,y) - h\x,y)\ < M^y exp(Mn), 
where M = max{M2,M5} (see (4.15)). 

Proof. Assume ?y > (the case y < is proved similarly). 
IfO = 2/^<2/<yi, 

rv 

h{x,y)=x + / a{h{x,s))ds 
Jy'o 

and 

rv 

h\x,y) = h\x,yl) + {y-y^o)a{h\x,yl^))=x+ / a{h\x,yl))ds. 

Then, using that tr is a Lipschitz function, 

rv 

\h{x,y)-h\x,y)\< / \a{h{x,s)) - a{h\x,ylj))\ds 
Jo 

<M \h{x,s)-h^{x,s)\ds + M \h^{x,s) - h\x,yl)\ds. 
Jo Jo 

(4.27) 

Since, for y^ < s < y{, 

\h'{x,s)-h\x,yl,)\ < \h'{x,yl,) + is~y'o)aih\x,yl))-h'{x,y'o)\ 
< |s-y^||a(/i'(x,y^,))| <M(yi-y^,), 

then 

M r \h\x,s) - h\x,y'o)\ds < M^y[ - ylf = M^rf. (4.28) 

By (4.27), (4.28), and Gronwall's lemma, for = yl < y < y[, 

\hix, y) - h\x, y)\ < M^f cxp{M{y - y^)) < M^rf exp(M(t/i - y^)) 

<M^rfexp{Mn). (4.29) 

We will prove by induction that for /c = 0, • • • ,'nl — 1, if yl. < y < 2/[_,_i, then 

\h{x, y) - {x, y) I < M Vf [exp(M(t/[+i - y^)) + . . . + exp(M(y[+i - y{))\ . 

(4.30) 

For A; = we have the result by (4.29). For y\<y < J/fc+i, 

rvi rv 
h{x,y) = x+ a{h{x,s))ds+ / a{h{x,s))ds 
Jy'o Jvi 

rv 

= Kx,yk)+ I (T{h{x,s))ds, 
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and 

rv 

h\x,v)=h\x,yi) + {y-yi)<j{h\x,yi))=h\x,yi)+ / a{h\x,yi))ds. 
By induction on fc, (4.30), and a being Lipschitz, 
\h{x,y) - h\x,y)\ < \Hx,yi) - h\x,yi)\ + f \<j{h{x,s)) - aih' {x,yi))\ds 
<M\f [cxp(M(2/i - y'o)) + ■■■+ exp(M(2/i - yi_,))] 
+M i r \h{x,s)-h^{x,s)\ds+ \h^{x,s) - h'{x,yi)\ds] . (4.31) 

Since, for yi < s < yi^-^, 

\h\x,s)-h\x,y[)\ = \h\x,yi) + {s~yi)a{h\x,y[))-h}{x,y[)\ 
<\s-yi\\a{h\x,yi))\<M{yi^,-yi), 

then 

\h\x, s) - h\x,yl)\ ds < M{yi+, - yif = Mrf. (4.32) 

•'y'k 

By (4.31), (4.32), and Gronwall's lemma, 
\h{x,y)-h^{x,y)\ 

<M\f [exp(M(y[ - y[)) + • • • + exp(M(y[ - y^i)) + l] exp {M{yk+i ' Vk)) 
=M'rf [exp(M(t/[+, - y^)) + • • • + exp(M(y[+i - yi))] . 

Then for all {x, y) E [— n, n] x [—n, n] and I > n, there is some k G {0, - ■ ■ ,nl — l} 
such that y\.<y< yk+n and by (4.30), 

\h{x,y) - h\x,y)\ < M\f [cxp(M(?yt+i - y'^)) + ■■■ + cxp(M(y^+i - yi))] 
< M^rf (k + l) cxp(Mn) < M^rfnl exp(Mn) 
= M^jexp{Mn). 

□ 

Next wc do the approximation of Y. 
We denote 

a„ = n-i/2+'^+^'(logn)^/2, (4.33) 
with P,S such that \H - 1/2| < /3 < 1/2, 0<(5</3and/3 + (5< 1/2. 

Proposition 4.8. Let Y and Y" be the processes given by (4-4) O'lT-d (4-5), res- 
pectively. Then 

pflimsup{||F-y"|U>a„}) =0, 

\ n— >-oo / 

where an is defined by (4-S3). 
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Proof. 



\Yt-Y,-\ 



dh 
\dxi 



(y„Bf) h{h{Y,,Bf))ds 



< [ Ii{s)ds+ [ l2{s)ds 
Jo Jo 



where 

his) - 
and 

his) = 



dh 



dh 



iyr,B^) 



As h is bounded by Mi, then 



/ his)ds <Mi[ 
Jq Jo 



dh 



dxi 



dh 



dxi 



Denoting 



Fixi,X2) = ( ^(^I'^z^ 



exp I — / a'ihixi,u))du 



(see (2.12)), we have 



— (Y B") 

dx^ ' 



dh 

— fF" 
dxi^^ 



< his) + his), 



where 



his) = \FiY,,B^)-FiYr,B^)\, 
his) = \FiY:\B^)-FiY:\B^)\. 

By (4.20), 

his) < M3||B^||ooexp(2M2||B^||oo)|n 
Since a' is bounded by M2, and by (4.17), 



dF 



exp 



a {h{xi,u))du 



(4.34) 



\bihiY:,B^))\, (4.35) 



|6(/i(y„ Bf )) - bihiYp, B^))\ . (4.36) 



ds. (4.37) 



(4.38) 



(4.39) 



|(t'(/i(xi,X2))| < M2exp{M2\x2\). 



By (2.8) there is measurable subset A of the underlying sample space with 
P(A) = 1, and for each 00 & A there is a positive integer N — Niui) such that 



|B^(a;) -B"(w)||oo < 1 for ah n>N. 



(4.40) 
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For each oj ^ A fixed and n > N{oj), by the mean value theorem, for some r(w) 
between (w) and -B"(a;), |f(w)| < 1 + ||B^(a;)||oo, thus (omitting w) 



his) 



dF 

— (Y" f] 



\Bf S;' I < Af2 cxp(M2|r-|)|i3f - 



< M2exp(M2(l + ||B^||oo))||S^-B"||oo for n>N. (4.41) 
From (4.37), (4.38), (4.39) and (4.41), for n>N, 

fh{s)ds < Ml AM3||B^|Uexp(2M2||B^|U)|n 
Jo Jo 

+M2exp(M2(l + ||B^|U))|-Bf - B^\]ds. (4.42) 

Now, by (4.17) and (4.36), 

his) < eMM2\\B"\Uihis) + his)), (4.43) 

where 

his) = |fe(/i(n,Bf )) - bihiY:,B^))\ , 

his) = |6(Mn",5f )) - bihiYr,B:))\ . 

Prom (4.21) and (4.22), 

his) < M4exp(M2||B^||oo)|n (4-44) 

76 is) < M4M5 |Bf -B^\. (4.45) 
By (4.43), (4.44) and (4.45), 

/ his)ds< / exp(M2||S«||oo) [M4exp(M2||B^||oo)|n-n"l 
Jo Jo 

+M4M5|Sf -B^ljds. (4.46) 

Therefore, for w G A there is iV = 7V(a;), as above, such that for n > TV, by (4.34), 
(4.42) and (4.46), 

|it-i"i"i < alls'" -s"iioo+« / in-n"Ms> 

Jo 

where a and k are the random variables 

a = (TMiM2exp(M2) +TM4M5)exp(M2||B^||oo), 

K = (MiMallB^lloo + M4)exp(2M2||B^||oo). 
By Gronwall's lemma 

\Yt-Yt^\<Z\\B" -B^W^, 
where Z is the random variable 

Z = aexpiTK). 

Then, for n> N, 

\\Y-Y-\\^<Z\\B"-B-\U 

and it follows that on A, 

limsup{||y - y"||<^ > C hmsup {Z\\B" - > "n} ■ 
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Hence 

P flimsup{||y - > an}) < P (limsup {Z\\B" - B^\\^ > a^}] • 

\ n—^oo / \ n— >-oo / 

(4.47) 

Next we show that 

hmsup{Z||B^ -B"||oo C limsup (llB^ -B"||oo > n-^/^+''(logn)^/H . 

Since 

hmsup {Z\\B" - > n-V2+^+^(logn)5/H 

= limsup {{Zn-^)\\B" - B"||oo > n-^/^+^ {lognf^^} , 

n— >-oo ^ 

then, by (2.8) we have 

P flinisup{Z||B^ - B"||oo > 

\ n— ^oo / 

The proof is finished by (4.47) and (4.48). □ 

Now we do the approximation of Y". 
Lemma 4.9. For every n,m and y"-"* defined by (4-9). 

\Yl'n < \xo\+TM,eMM2 ||P"IL), t G [0,T]. (4.49) 
Proof. As b and a are bounded, and by (4.17), then, for /" defined by (4.8), 

|r(y,"''",Br)| < Miexp(M2 mj. (4.50) 
We will first prove that for all fc = 0, • • • , m — 1, 

< \xo\+krmM,eMM2 \mU. (4.51) 
For fc = it is obvious. If we assume that 

IFttTI < ko| + (fc-l)r„Miexp(M2||P"||J, 
then from (4.50) and (4.9), 

< \xo\ + {k- l)r„Miexp(M2 ||P"||^) + r„Mi exp(M2 ||S"||^) 
= |a;o| + fcr„Miexp(M2||B"||^). 
Now, iitk<t< tk+i, by (4.50) and (4.51), 

= \Yrr+it-tk)nYrr,Bi)\ 

< \xo\ + kr„,M, cxp(M2 ||P"IL) + r„Mi exp(M2 ||B"|U) 
= \xo\ + (fc + l)r„Mi exp(M2 ||P"||^) 

< |a;o|+mr„Miexp(M2||P"||^), 

and (4.49) is obtained. □ 
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Lemma 4.10. There exists N = N{uj) > such that for every n, m and t G [0, T], 
(yn,m^ B,") G [-N, N] X [-N, N] a.s. (4.52) 

Proof. Let w e ^ and N = N{u}) (see (4.40)), then for all n> N, 

||S"||oo < l + ||B^||oo (4.53) 
and by the Lemma 4.2, for n < N, 

\\B"\\^<KN^+^T. (4.54) 

Hence, by (4.49), 

lyn.m, < [ko|+TMiexp(M2(l + ||B^||oo)) ifn>7V, 

'* ^ - \\xo\+TMiexp{M2KN^+^T) if n < TV. ^' ' 

Taking 

iV - iV (a;) - max | ^^^^ ^ exp(M2i^7Vi+/5T) / ' 

(4.56) 

we then have the result by (4.53), (4.54), (4.55) and (4.56). □ 

Corollciry 4.11. Let N be given by Lemma 4. 10. Then for all n > N and t G 
[0,T], 

\h{yn,m^ ) - /i"(y,"''", B^)\ < M^^ exp {MN), 
where M = max {M2, M5}. 

Proof. Follows from Lemma 4.7 and (4.52). □ 
Lemma 4.12. For every n > N, m and t G [0,T], 

\Yr - F"'™! < Jn,mT (exp (ZiT)) a.s., (4.57) 
where Jn,m is the random variable 

Jn,m =ZiM, exp(M2 ||B"||^)r„ + Z2Kn^+^rm 

+ M2exp(M2||B"||oo)exp(M7V)(MiM3||B"||oo + M4)-, (4.58) 

n 

and Zi,Z2 and N are the random variables defined by (4.25), (4-26) and (4.56) 
respectively, and rm = T/m. 

Proof Firstly, if = io < t < ti, by (4.5), (4.7) and (4.9), 

IV - Yr^\ < f \f{Y-,B-) - f-{Y;i^^,Bl)\ds, (4.59) 

J to 

and 



i/(y,",i?:')-r(r,:'",i?r„; 



<i/(n",K) - /(n"'",5")i + i/(n"'",K) - f{YtT^B^)\ 

+ \f{Yrr^B:) - f{Y,Y\Bl)\ + \f{Y-^"\Bl) - f-{Y--"',Bl 
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By part 2 of Lemma 4.6, 
Analogously, and by (4.9) and (4.50), for s <ti, 



(4.61) 



< Z,{s-to)\nYr''^,Bl)\ < Z,M,exp{M2mjrm. (4.62) 



'to J -^s l\ — ■^i-l-^ s -'to 
"^to '-°toJ 

By Lemmas 4.2 and 4.6, for to < s < t\, 

|/(y,^'™, - fiYtT^ K)\ < Z^\B: - BZ I < Z2Kn'+^{s - to) 

< Z2Kn^+^rm. (4.63) 

Moreover, 

l/(^t:'™, BZ) - nVr, BZ)\ <Ir + h, (4.64) 

where 



Ir = \b{h{Yrr,BZ))\ 



exp 



a'(/i(y"''",n))rf« 



— exp 



a'{h"{Y,"^'"',u))du 



/« = 



exp - / c7'(/i"(y,;:'™,u))du 



\bih{Y,"'"^,Bl)) - 6(/i"(y"'™,B"))l- 



Applying the mean value theorem for the exponential function and a' , by Lem- 
mas 4.7 and 4.10, (4.17), and boundedness of b, a', then, for n> N, 



Ir < MiM3exp(M2||B'^||oo) 



\h{Y-'"^,u)-h-{Y-'"^,u)\du 



< MiM3exp(M2||B"||oo) / — exp {M N)du 

Jo ^ 

< MiM3M2exp(M2||B"||oo)exp(M7V)||B"|U-. 

n 

By Corollary 4.11, (4.17) and b being Lipschitz, for n > N, 

h < M4exp(M2||B"||oo)|MC"'St")-/^"(^t:''",St"o)| 

- N 

< M4M2exp(M2||B"||oo)exp(MAr) — . 

n 

Prom (4.64), (4.65) and (4.66), 

\f{Y,Y'\BZ)-nY,Y\BZ)\ 

< M^exp{M2\\B"\\^)cxpiMN){MiM3\\B^^ + Mi) — . (4.67) 

n 

By (4.60), (4.61), (4.62), (4.63), (4.67) and (4.58), for s £ [0,T], 

y-n,m r>n 
■^to ' -"to 



(4.65) 



(4.66) 



|/(F,", B^) - r(C™, Bl)\ < Z,\Y^ - y,"'-| + J„,„. 



(4.68) 



20 



JOHANNA GARZ6n, LUIS G. GOROSTIZA, AND JORGE A. LE6n 



From (4.59), (4.68), for < t < h, 

J to J to 

and by Gronwall's lemma, for to <t <ti and = T/m, 

\Yr - F"'™! < Jn,mrm exp (Zi(ti - to)) . (4.69) 

If ii < i < t2, 



Jti 

< lYt: - Yt^n + f \f{Y:, B2) - r{YtT^ K)\ds- (4.70) 

Proceeding similarly, taking ti instead of /;o, from (4.68), 

\f{Y:,B-)-r{Y^T,B-^)\<Z,\Y:-Y^^"^\ + Jr,,m. (4.71) 
From (4.69), (4.70) and (4.71), for ti<t< t-2, 



Jn,mrmexp{Zi{tl-to))+ [ Zi - F,"''" + / Jn,mds 

Jtl Jtl 



< Jn,mrm{exp{Zi{h-to)) + l)+ - F,"'™ |ds. 



Jtl 



and by Gronwall's lemma, 

< J„,„r„ (exp (Zi(ti-io)) + l)cxp(Zi(i2-ti)) 

< Jn,mrm (exp (Zi (t2 ^ ^o)) + exp {Zi{t2 - ti))) 

< Jn 2(exp(ZiT)). 
Analogously, for tk < t < tfc+i, fc = 0, ■ • • , m — 1, 

<Jn,mrm[exV (^i(tfc+i - to)) + exp {Zx{tk+i - ti)) H h exp {Zi{tk+i - tk))] 

<Jn (fc + l)(exp(ZiT)) < (exp(ZiT)) 
=Jn,mT (exp(ZiT)) , 
which finishes the proof. □ 

Proposition 4.13. Let and y"'"* be given by (4-5) and (4-9), respectively. 
Then 

P(limsup|||y"-y"'"'||oo>a«}) =0, 
where an is given by (4-33). 

Proof. By Lemma 4.10, there is A'' = ^(0;) given by (4.56) such as for n, m, 

(y,"'", B^) G [-N, N] X [-N, N] a.s. 

The random variables Zi and Z2 given by (4.25) and (4.26) are bounded uni- 
formly in n, and defining M = max {Mi, M2, M3, M4, M5}, 

Zi < exp(2MA)[M2iV + M], Z2 < 2M^ exp(MiV), 
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and by (4.58), 



Jn,m <cxp{3MN){M^N + M^)rm + 2M'^Kexp{MN)n'^ 



N 



+ ]Vrcxp{2MN){M^N + M) — . (4.72) 

n 

From (4.57) and (4.72), for n > TV and taking m = n^, then rm = T/n^ and 

|>;"-i;"'"'|< Jn,mr(cxp(ZiT)) 

< {exp(3MAf)(M3Ar + M'^y^ + 2M'^K exp{MN)n^+'^rm 

exp{2MN){M'^N + M)^^Texp {T exp{2M N)[M^ N + M]) 



< Z3- ^ 



where the random variable is defined by 
Z3 = [ex:p{iMN){M^N + M'^)T + 2A'PKkxp{MN)T 

+M'^N exp(2MiV)(M^ at + M)}T exp (T exp(2MAf ) [M^AT + M]) . 
Therefore we obtain 

P(hmsup{||y" - F"'"'||oo > an}) 

<p(limsup{Z3 > n^-f^an}) = p(limsup{Z3 > n^^^+\\ogn)^/'^}) = 0. 

□ 

Corollary 4.14. LetY andY"'"^ be given by (2.11) and (4-9), respectively. Then 
p(limsup{||y-y"'"'||oo >a„}) =0, 

where an is given by (4.33). 

Proof. By the proofs of Propositions 4.8 and 4.13, and replacing Z and Zs by 2Z 

and 2Z3, respectively, 

p(limsup{||F-y"'"'||oo >Q„}) 
<p(limsup{||y-y"|U >a„/2}) +p(limsup{||y"-y"'"'|U >a„/2}) =0. 

□ 

We define for each n = 1, 2, . . . , 

Xr = M^",Sn, (4.73) 
where P" and Y" are given by (2.6)-(2.7) and (4.5), respectively. 

Proposition 4.15. For any C > 0, 

p(limsup{||X - X"||oo > Can}) = 0, 

where X is given by (2.9) and an is given by (4-33). 



22 JOHANNA GARZ6n, LUIS G. GOROSTIZA, AND JORGE A. LE6n 

Proof. For convenience of notation we put C — I. Then 

\Xt - Iri = \KYt, ) - h{Y,-, B^)\ < I^{t) + ho{t), (4.74) 

where 

h{t) = \h{Yt,B^)-h{Yr,B^)\ 
ho{t) = \h{Y,-,B^)-h{Yr,B^)\. 

By (4.19), 

hit)<Z4\Y-Y^^, (4.75) 

where = exp(M2 ||B^||^), 

Proceeding similary as in the proof of (4.48), then from Proposition 4.8, 

p(hmsup{Z4 ||y - > an/2}) 

< P(hmsup{||y-F"||oo >n-l/2+/3+5/2^jQg^^5/2|\ 

Thus, by (4.75), 

P(limsup{||/9||cx, > an/2}) = 0. (4.76) 
Prom (2.10), boundedness of a, and the mean value theorem, 

/io(i) <M5||B^-P"|U, 

and by (2.8), 

p(limsup{||7io|U > an/2}) = 0. (4.77) 

The result follows from (4.74), (4.76) and (4.77). □ 

Por the final step of the proof of the theorem we go to (4.11). 
By (2.9), (4.73), (4.12) and Proposition 4.15, 

p(limsup{||i?i||oo > an/4}) = 0. (4.78) 

By (4.52) there is N = N{uj) given by (4.56) such that for each n, m, 

iY,"-"',B^)€[-N,N]x[-N,N] a.s. (4.79) 
Prom (4.13) and (4.19), 

H2{t) < exp(M2||i?"||oo)|i^t" - < Z^lYr - (4.80) 

where Z5 = exp{M2N). By Proposition 4.13 and (4.80), 

p(limsup{||if2||oo > an/4}) < p(Hmsup{4Z5||F" - F"'"' ||oo > a„}) 

< P(limsup{||y" - y"'"'||cx> > n-^/^+^+^/^{lognf/^}) = 0. (4.81) 

^ n— >oo ^ 

On account of (4.79) and (4.14), for all n > N, applying Lemma 4.7, 

Hsit) < Ze-, 
n 
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where Ze = PPNcxp{MN). Then, 
p(limsup{||iJ3||^ > an}) < PUimsupIZe > na, 

< p(hmsup{Z6 > nV2+^+^(logn)5/2}) = 0. (4.82) 

The final result follows from (4.11), (4.12), (4.13), (4.14), (4.78), (4.81) and 
(4.82). 
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